Magnetic collapse of a neutron gas: Can magnetars indeed be formed? 
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A relativistic degenerate neutron gas in equilibrium with a background of electrons and protons 
in a magnetic field exerts its pressure anisotropically, having a smaller value perpendicular than 
along the magnetic field. For critical fields the magnetic pressure may produce the vanishing of the 
equatorial pressure of the neutron gas. Taking it as a model for neutron stars, the outcome could 
be a transverse collapse of the star. This fixes a limit to the fields to be observable in stable neutron 
star pulsars as a function of their density. The final structure left over after the implosion might 
be a mixed phase of nucleons and meson condensate, a strange star, or a highly distorted black 
hole or black "cigar", but no any magnetar, if viewed as a super strongly magnetized neutron star. 
However, we do not exclude the possibility of a superstrong magnetic fields arising in supernova 
explosions which lead directly to strange stars. In other words, if any magnetars exist, they cannot 
be neutron stars. 

INTRODUCTION 

Gravitational collapse occurs in a body of mass M and radius R when its rest energy is of the same order of its 
gravitational energy, i. e., Mc 2 ~ GM 2 /R. We would like to argue that for a macroscopic magnetized body, e.g., 
composed by neutrons in an external field B {\B\ = B), new physics arises and a sort of collapse occurs when its 
internal energy density U is of the same order than its magnetic energy density M. ■ B, where M. is the magnetization. 
This problem is interesting in the context of both cosmology and astrophysics, as for instance in the study of objects 
such as neutron stars. A gas of neutral particles having an anomalous magnetic moment (as a model for neutron 
stars (NSs); here we assume, as usually, a background of electrons and protons in j3 equilibrium, which is demanded 
by Pauli's Principle to guarantee neutron stability), when placed in extremely strong magnetic fields has a nonlinear 
(ferromagnetic) response to the external field and is also unstable due to the vanishing of the transverse pressure for 
surface fields strong enough (B sur f > 10 16 G). In this phenomenon quantum effects play an essential role due to the 
coupling of the particles' spins to the microscopic field B seen by the particles (spin-polarization). For fields of this 
order of magnitude there are values of the density for which the magnetic energy of the system becomes of the same 
order of magnitude than its total energy. At these physical conditions any structure of super dense matter composed 
of neutral particles having a magnetic moment may undergo a transverse collapse when its pressure perpendicular to 
B vanishes. This implosion is driven by the same mechanism described in Jl| for charged particles. 

We present in this new paper, which is a more elaborated version of |2|, the main ideas concerning the role of 
ultra strong magnetic fields in a gas of neutral particles: The standard model of NSs, as envisioned by Duncan and 
Thompson in their model of magnetars^. The fundamental result that we obtain shows that a NS, i. e., a neutron 
gas permeated by a super strong magnetic field is unstable and must collapse. This inedit result seems to ban the 
possibility of formation magnetars. We stress in this respect the fact that similar results were obtained by two different 
groups: Khalilov 4] and Ghosh, Mandal and Chakrabarty^. Although both teams of researchers arrived to the same 
conclusion like ours, quite different approaches were pursued. 

The paper is organized as follows: Section II reviews the concept of anisotropic pressure in self-gravitating systems 
like NSs. In Section III are given the tools to construct the energy- momentum tensor of a neutron gas. Section IV 
discusses the main differences between the classical and quantum collapse of a gas configuration in an approximation- 
independent way, based on the sign of the electromagnetic response of the medium to the external applied field. In 
Section V the dynamics of the neutron gas composing a NS is studied in the presence of a magnetic field, followed in 
Section VI by the derivation of its thermodynamical potential and magnetization. Section VII discusses the conditions 
for the collapse to take place, and an application of this physics to the stability analysis of the proposed magnetars 
is given. The main result of this paper shows that such objects should not form if they are envisioned as a standard 



2 



neutron gas in Fermi beta equilibrium as is claimed in the original idea introducing the concept of magnetar. Some 
closing remarks and further potential applications of this theory are part of the final section. 

ANISOTROPIC PRESSURES IN SELF-GRAVITATING SYSTEMS 

The issue of local anisotropy in pressures was extensively reviewed by Herrera and Santos 6] in a general relativistic 
approach. These authors present several physical mechanisms for its origin in both extremely low and very high 
density systems, which may include astrophysical compact objects. In the case of highly dense systems, it was 
pointed out that "exotic" phase transitions could occur during gravitational collapse, which is exactly the problem we 
are concerned with in this paper. However, the more fundamental problem regarding the appearance of anisotropic 
pressures in strongly magnetized compact stars was left open. More recently, Mak and Harko Q present a class of 
exact solutions of Einstein's equations corresponding to anisotropic stellar configurations which can describe realistic 
neutron stars. 

We want to provide a more detailed description of the arising of anisotropic pressures in a relativistic system like a 
NS, an essential point in understanding the physics behind the problem of stability of ultra magnetized compact stars. 
We shall give firstly general arguments to support this our view, and then we concentrate in the specific calculations, in 
the one-loop approximation, of the thermodynamic potential of the neutron star configuration and its magnetization, 
the properties upon which the most crucial conclusions can be drawn. 

To fix ideas, we shall work in the grand canonical ensemble, and we are considering as subsystem, some region 
inside the star. Such subsystem is under the influence of the magnetic field created by the rest of the system, which 
we name H (\H\ = H). The response of the subsystem is to polarize itself creating a magnetization A4 in the medium 
(the neutron gas) satisfying the relation: H = B — 4ttA4. Obviously, inside the subsystem the microscopic field is 
B = H + 4ttA4, since B (named also magnetic induction) and the electric field E are the true fields actingon the 
electric charges and magnetic dipoles of elementary particles [|| 0- (Note, however, that Landau in Ref. 10] uses 
the notation H to name the magnetic field in vacuum). The field B, as well as E, satisfy the Maxwell equations for 
particles in vacuum. In what follows, when we are to discuss the dynamics of the particles in the neutron gas we will 
sometimes refer to B as the external magnetic field, as is usually named in quantum field theory and astrophysics. For 
an external distant observer, B = H (in Gaussian units) since the magnetization is assumed to exist only inside the 
star. We emphasize that actually B and H are external fields within different contexts: H is external to the subsystem 
object of study in the grand canonical ensemble, whereas B is external to any particle chosen in the subsystem (it feels, 
in addition to H, the contribution from the magnetization field AirM due to the other particles of the subsystem). 

In the case of a gas of electrically charged particles in an external constant magnetic field B, in classical electro- 
dynamics, it is the Lorentz force F = ev x B/c the source of an asymmetry in the pressure components parallel and 
perpendicular to B. By writing ev = jAV, where Ay = dxidx2dx%, calling /j = Fi/AV as the i-th component of 
the force density, and substituting j = cV x M, one has 

fi = ~(d l M s )B s + {d s Mi)B s . (1) 

By multiplying by AV^ = dxydx^dx^ and assuming B s = BS S 3 and dAii/dx3 — (actually it is also M.i — AiS^), 
only the first term in Q remains nonzero and one recovers an expression for the Lorentz force, which is obviously 
perpendicular to the field B. For the corresponding pressure it yields Pl± = —M ■ B. This is a classical effect and 
obviously Pli_ must be added to the usual kinetic isotropic pressure Pq, so that the total transverse pressure becomes 
P± = Pq + Pl±- As in classical electrodynamics, by Lenz law, M. is opposite to B (spin effects are neglected), then 
M. ■ B < 0, and Pl± > 0. The opposite case occurs when M. is parallel to B, which occurs in the quantum case, i.e., 
when spin effects are taken into account. We alert, however, that in the definitions and derivations that follows use 
will be done of the magnitudes \B\ and \M \ = A4 of both vectors B and M. instead of the vectors themselves. 

THE ENERGY-MOMENTUM TENSOR OF A NEUTRON GAS 

Based on more fundamental grounds, one may write the general structure of the energy-momentum tensor of a 
neutron gas in an external field B in the same way as one does to construct general tensors, as the polarization 
operator tensors [ill ], for instance. In an external field F^, in addition to the basic 4- velocity vectors of the medium, 
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u^, and particle momentum we have two extra vectors F^vky, F 2 k v , to form a basis of independent vectors (in 
what follows we shall use the notation = F^\F\ V ) . From them we may build a set of basic tensors which, together 
with the tensors 6^, F^, F^ v , serves as a basis in terms of which we can expand any tensor structure related to the 
particle dynamics, in particular, the energy-momentum tensor. However, to get rid of tensor structures containing 
off-diagonal terms, which would correspond to unwanted shearing stresses in the fluid rest frame u — (0,0, 0,1*4), we 
exclude some of them, i. e., /c M fc„, k v F^\k\, F^ v , u^k v , or any of its combinations. By following the arguments used in 
we conclude that we are left in the present case with three basic tensors: 6^, F^ u , u^u v to describe the dynamics 
of such a neutron gas. Thus, the structure of such an energy-momentum tensor is then expected to be of the form 

= adfj, v + bF 2 v + cu^Uy (2) 

where /1, v = 0, 1, 2, 3; and a — p is the isotropic pressure term, b = M./B and c = U + p. In the present case the 
second of these tensors can be written in a simpler form as F 2 U = —B 2 8j iv . The tensor (J2J has then the spatial 
eigenvalues p — BM.,p — BM,p 1 and the time eigenvalue — c = —U — p, since u^u^ = — 1. These eigenvalues exhibit 
the anisotropy in pressures perpendicular and parallel to B. 

In dealing with a quantum gas in an external field we shall assume that the sources of the field H are either 
classical (currents) or due to quantum effects. Although it is out of the scope of this paper to discuss the mechanism 
for producing such field, we suggest, however, some viable sources able to induce a self-consistent field, as for instance 
a condensate of the vector meson p, the neutron spin-spin ferromagnetic [l^ coupling (see below) , or even by diquark 
[l3j | condensation. 

As in the case of the electron in an atom the basic dynamics in our present case is described by the Dirac equation 
in an external field, in place of the Lorentz force. According to the Ehrenfest theorem, classical dynamics is contained 
on the adequate average of quantum dynamics. But quantum dynamics leads also to several new phenomena not 
having classical partner. After solving the Dirac equation one gets the energy eigenvalue spectrum |l4|. This energy 
depends on the microscopic magnetic field B through some interaction term in the initial Lagrangian. These energy 
eigenvalues, after the quantum statistical average, determine the thermodynamic properties, such as the neutron gas 
pressure. If the coupling constant is turned to zero, the particles would not feel any pressure coming from the external 
field. Thus, since the classical Maxwell stress tensor of the field H does not express by itself the interaction with the 
particles, i. e., it expresses the momentum and energy of the external field, we have no need to add it to the expression 
(@J below. The tensor Q contains already the basic tensor structures of the problem, including the Maxwell stress 
tensor of the field B, which depends on the external field H and the magnetization \ A4\. 

The total external field B contributes with virtual particles, expressed through the Euler-Heisenberg vacuum terms 
arising in the regularization of the quantum vacuum terms |15j . These vacuum terms appear in the calculation of 
the basic statistical quantity, the thermodynamical potential: fl = — (3~ 1 lni?, where the grand partition function, 
is built up on the particle spectrum. Our thermodynamical potential is the sum of two terms, f2 = fl s t + flv, the 
finite statistical term fi st plus the vacuum field contribution fly, which is divergent. In the process of regularization, 
it absorbs the classical field energy density term B 2 /8tt. We observe here that Landau, in p. 69 of Ref. 16], uses the 
specific term thermodynamic potential, denoted by $ in that reference, to name what in western literature is known 
as Gibbs free energy, which is denoted by G. Our thermodynamic potential is just what Landau defines as "new 
thermodynamic potential CI" , but we have taken it per unit volume, which is dependent on T and /1 in absence of 
external field, and in our case is dependent also on B. Observe that fl = F — G. In the zero field case, it would be 
n = —P, where P is the isotropic pressure. Due to the spatial anisotropy introduced by the magnetic field B, the 
pressures are not the same in all directions, and only in the direction parallel to B it adquires the value Cl(B) = —P3 
(See below). 

The coupling of the spin dipole moment of neutrons in an external magnetic field B produces a loss of rotational 
symmetry of the particle spectrum (in what follows we will consider B along the X3 axis). From the spectrum, 
which is expressed in terms of B, by the standard methods of finite temperature quantum field theory, we obtain the 
thermodynamical potential (per unit volume) CI = Cl(B) 1 and from it all the thermodynamic properties of the system, 
in particular its magnetization, as is done by Huang [ljj (p. 237), which is the statistical average M = — (dCl/dB). 

Notice that this definition is consistent with our convention by which thermodynamical quantities are defined in 
terms of the microscopic magnetic field B acting on the particles as an independent variable (see for instance Ref.||), 
instead of using the quantity H. Thus, the thermodynamical variable conjugated to the magnetic field B is the 
quantity A4(B), the system magnetization as introduced above. We may write then 



Cl = - MdB - P , 



(3) 
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where Pq = — 0(0) is the term corresponding to the zero magnetic held pressure. One must emphasize that in the 
quantum relativistic case Q depends on B nonlinearly. From the explicit expression for £1 given in another section 
below, one finds that the dependence of the energy spectrum on the particle momentum is not rotational invariant. 
This fact determines a reduction of the symmetry of the otherwise isotropic thermodynamic properties of the system 
such as the pressure, which is expected to be axially-symmetric for the reasons pointed out above. 

By using the Green functions method it is found that the energy-momentum tensor of matter in an external constant 
magnetic field obeys the general structure J5J), and we have 



V = (Tdn/dT + J2^rdn/d^ r )s 4 ^ (4) 

+ iF^Fuxdfl/dF 2 - s^n, 

where r run over the species involved. Below we will take r = n,p, e to describe teh neutron, proton and electron 
component of the star gas. Expression Q in the zero field limit reproduces the usual isotropic energy-momentum 
tensor = P5 M „ — (P + U)S4^5 V 4 of a perfect fluid. From Eq.@ the spatial components are 733 = P3 = — fi, T\\ — 
T22 = P±_ = — — BM.. The time component 7oo = — U r = —TS r — u r N r — il r , in which U r ~ fi r N r ~ Q r are 
quantities of the same order of magnitude. We shall assume that inside a stable NS there is locally a balance between 
the mechanical stresses (and forces) from the pressure exerted by the quantum gas and the gravitational stress, i. e., 
the gravitational force per unit area, exerted by the star mass, as defined below by equations ©,10. By abusing a 
bit of the phrasing we shall refer over the paper to this last one as the "gravitational pressure" , in an attempt to turn 
more meaningful discussions regarding effects of the gas pressure and those stemming from its own gravity. In the 
case of a neutron star the correct treatment of this problem requires of general relativistic corrections 18J. However, 
for the general purpose of this paper we consider here the problem in a rough manner, in a similar fashion as it was 
done in the white dwarf case with zero magnetic field , in Newtonian gravity. If dSi is an element of surface and 
dxi an increase in the coordinate z = 1,2, 3, one can express the balance between the work done on the star and the 
variation of its gravitational energy, E g = GM 2 R~ X , under a small change of volume AV as 



/rt W 
(-TijdSj + —^)d Xl = , (5) 
OXi 

where the quantity in parenthesis must be zero, and describes the balance of forces. We will use now cylindrical 
coordinates with the z axis parallel to B, and assume that the shape of the deformed star can be (very roughly) 
approximated by the simplified model of a cylinder of height Z and radius r±_ (we could equally use as a simple 
model a cylinder of height Z ended by two half-spheres of radius r± . The equations © and J7J below would differ 
by unimportant numerical factors). We assume also that the change of the cylinder's volume AV and surface AS are 
small, so that the total surface of the star 5* = 2irr±(r± + Z) is approximately constant. Then we have dS±_ = 2?rrj_dz, 
dSs = 2nr±dr± and integrating we get the equilibrium between the gravitational and gas pressures P±_ = P g ± and 
= Pg33, where 



and 



9± - ~ 2nr ± Z dr ± ' W 



P 1 9E ° (7) 

^33-—^- (7) 



ur 



Eqs.lO,® are to be interpreted respectively as the transverse and longitudinal "gravitational pressures" in the sense 
defined above. By assuming that -g— 1 and -^f- are quantities of the same order, to preserve the balance of transverse 
pressures when P± decreases (due to the increase of the product BAi), the quantity l/2Trr±Z must decrease. This 
is achieved whenever Z increases faster than r± decreases (all this can be seen by assuming that the total surface 
S = const and dr±/dr < 0, where r is a parameter like time). In that case the outcome is a stretching of the body 
along the direction of B. Thus, the anisotropy of pressures in our problem suggests that matter in the body is also 
distributed anisotropically, what leads to prolate isobaric surfaces j^Qj . 

The condition P± = implies also that U ~ BAi. As P3 > 0, an instability arises in the system leading to a 
transverse collapse. Thus, as pointed out in the introduction, new physics arise in connection with the vanishing of 
the transverse pressure, P± = 0. This peculiar behavior will be discussed in the forthcoming sections. 
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By writing A4 = (B — H)/4ar, one may formally write il = — -^B 2 + J HdB — Pq. We remind that as 
il = F — G, the last expression is consistent with what would be obtained in the classical nonrelativistic case , 
where F = F + j HdB/4ir is the Helmholtz free energy and G = F + j MdB + P Q = G + B 2 /8ir as the Gibbs free 
energy. Because of our definition of M., the last term is given in terms of B and not in terms of H. As we have 

T 33 = -ft = j-B 2 - i- / HdB + P , (8) 

and also 



T n =T 22 = -ri-BM = -^-B 2 + ^- [ BdH + P , (9) 

it is straightforward to see that the spatial components of the energy-momentum tensor %j = 1,2,3) can be 
rewritten as 



%j = P % - %f (B, H) + S{B) l3 , (10) 

where S{B)ij = -^[BiBj — ^(B 2 )Sij] is the Maxwell stress tensor for the microscopic field B, and (B) = 
-^[HiBj — (J BdH)Sij] is the Minkowski tensor for relativistic nonlinear media (see below). This last term reduces to 
the usual expression in the nonrelativistic limit if H depends linearly onB 0. The definition (fLU|) expresses the total 
pressure as a sum of an isotropic pure mechanical pressure (independent of the electromagnetic field) plus a pressure 
coming from the Minkowski tensor due to the interaction of the external field H with the microscopic field B, plus 
the Maxwell stress tensor of the microscopic field B. If H = 0, B = 4ttA4 and %j = PoSij + S(B)ij. In this case the 
magnetic field is kept self-consistently. 

At this point we want to refer to the recent paper by Khalilov where a similar problem to the present one is 
studied. However, the expression for the stress tensor in a medium is taken as given only by the linear approximation 
of the Minkowski tensor term. This approach is not justified in the relativistic case. Further, the Maxwell tensor of 
the microscopic field B and the isotropic PoSij terms are omitted. This leads the author to wrongly conclude that the 
collapse occurs like in the classical case (see below), that is, along the B field, in contradiction to our present results 
and those of 0, 0- A consistent approach, as we have followed here, and discussed in the accompanying paper [20| . 
leads to opposite results compared to those of Ref. Q in what concerns the spatial direction of collapse, although the 
fundamental issue regarding the collapse of the neutron star remains taking place in both theories. 

Note, in addition, that Refs. 0and0, where the problem was also studied, did not take into account the dynamical 
effects of the strong (and super strong) magnetic fields supposed to exist in the core of canonical NSs, see for instance 
Refs. 0,H2, H3- H3- Thence, it is the contention of this paper to address this open issue. The novel results 
obtained here point out to the occurrence of new physics and processes in the (relativistic) astrophysics of compact 
objects that were not manifest in previous papers. 



CLASSICAL VS. QUANTUM COLLAPSES 

In Ref. 0] we found that a relativistic degenerate electron gas placed in a strong external magnetic field B is confined 
to a finite set of Landau quantum states. As the field is increased the maximum Landau quantum number is decreased 
favoring the arising of either a paramagnetic or a ferromagnetic response through a positive magnetization M. 1 up 
to the case in which only the ground state is occupied. The gas then becomes topologically one dimensional, and in 
consequence the pressure transverse to the field vanishes for fields B = VI/ ' M. 1]. Thus, the electron gas becomes 
unstable due to the decrease of the transverse pressure for fields strong enough, and the outcome is a collapse. 

For neutrons the magnetization is always positive (see arguments below) and nonlinear, what leads to a sort of 
ferromagnetic behavior. For fields strong enough the pressure transverse to the field, P±_ = — O — BM., is considerably 
decreased and may vanish. If we assume that extremely magnetized NSs, as the Duncan and Thompson magnetars|3J, 
have fields H ~ 10 15 G, and that inside the star B increases by following a dipole law B(r') = B sur f/r^, we expect 
near its surface magnetic fields ranging from 10 16 — 10 17 G|3j,|2j| up to values of order 10 20 G in its core|23j, where the 
field is maintained self-consistently, i. e., H = 0. For fields of this order of magnitude super dense matter composed 
of neutral particles having a magnetic moment may undergo a transverse collapse since P±_ vanishes. As discussed 



6 



below, the emerging physics seems to ban the possibility of magnetar formation. The outcome of such a collapse 
might be a compact star endowed with canonical magnetic field, as discussed below. 

In the classical case in which the response of the medium is due to the Lenz law, the magnetization is opposite to 
the external field H and it may happen that Ai < 0. This also occurs in the diamagnetic case. Then H > B and 
Pi > P3. Note that the opposite occurs in some permeable materials where Ai > and H = B — AnAA is small in 
comparison to either Ai and B; this is due to ferromagnetic effects which have quantum origin, as in the neutron gas. 

In the case of a classical magnetized gas, as P± > P3, this leads to the Earth- like oblatening effect described above. 
But opposite to this, for the critical quantum configuration of the NS gas the coupling of the spin magnetic dipolc 
with the magnetic field B plays the main role, and Ai > (see Eq.(j2U below and the subsequent discussion where 
this is shown explicitly), leading to ferromagnetic effects. The situation then is reversed and Pj_ is smaller than P3 in 
the amount BAi and it vanishes for -^B 2 — J BdH + Po leading, conversely, to a prolate configuration. 

In classical electrodynamics 01 ^ is suggested that the total pressure is given by the sum of the Maxwell stress 
tensor 5 M „ plus an isotropic pressure (Po) term. In the case of a constant magnetic field parallel to the X3 axis, the 
total pressure tensor reads %j = Po% + Sy or P 3 = P — B 2 /8n and P± = Pq + B 2 /8n. 

As pointed out before, the body deforms under the action of these anisotropic pressures. If the longitudinal pressure 
decreases, the body flattens along the magnetic field j^. Thus, in this pure classical case, for the extreme limit of 
flattening, P 3 = and P± = P + B 2 /8ir, and the body would collapse as a disk or a ring perpendicular to the field. 
Starting from general relativistic considerations it has been reported |22( the existence of a maximum magnetic field for 
having stationary configurations of NSs. (We interprete that result as indicating the occurrence of a classical collapse 
for fields larger than the maximum) . This field induces a toroidal configuration, which is topologically equivalent to a 
ring. In the quantum case, for degenerate fermions, as Ai > 0, it is Pj_ = — O — BAi which is decreased by increasing 
B. As the NS is in equilibrium under the balance of neutron and "gravitational" pressures, the body stretches along 
the direction of the magnetic field. Thus, for any density there are values of the field B high enough such that these 
pressures cannot counterbalance each other leading to a collapse perpendicular to the field for P± = 0. This collapse 
would leave as a remnant a nucleons plus a Bose-Einstein-like condensate, a hybrid or strange star with canonical 
magnetic field j^, or a distorted ("cigar-like") black hole. 

Our previous considerations are approximation-independent. In order to discuss an specific model, we shall start 
by computing the free particle spectrum for neutral particles in a magnetic field. 

THE NEUTRON GAS IN A MAGNETIC FIELD 

For free neutrons in a magnetic field B we have the Dirac equation for neutral particles with anomalous magnetic 
moment 

+ m + iqa^xFx^ip = 0, (11) 

where = |(7 m 7a —7a 7^) is the spin tensor, and P\ M is the electromagnetic field tensor describing B. By solving 
this equation we get the eigenvalues 



E n (p, B, rj) = \Jp 2 + (yfe + ml + m B) 2 , (12) 

where P3, p± are respectively the momentum components along and perpendicular to the magnetic field B, m n is 
the neutron mass, q = 1.91M„ (M n is the nuclear magneton), r\ = 1,-1 are the 03 eigenvalues corresponding to the 
two orientations of the magnetic moment (parallel and antiparallel) with regards to the field B. Here we make the 
assumption that the magnetic moment remains constant for large magnetic fields. Actually, this is not so. Radiative 
corrections change its value as a function of B. This problem is, however, beyond the scope of the present approach. 

From (|12J) we see that although the Hamiltonian is linear in B, the eigenvalues depend on B nonlinear ly. This 
makes the relativistic thermodynamic and electromagnetic properties of the system of neutrons very different from 
the nonrelativistic case. For instance, all thermodynamic quantities, Q, M. (and in consequence, H = B — 47rA / t), 
and N are also nonlinear functions of B. The expression (|12(l also shows manifestly the change of spherical to axial 
symmetry with regard to momentum components. This anisotropy in the dynamics is expected to be reflected in an 
anisotropy in the thermodynamic properties of the system, as it is expressed by the difference between the transverse 
and longitudinal pressures discussed earlier on an approximation-independent basis. 
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The partition function Z = Tr(p) is obtained from the density matrix describing the model p = 
e -Pf d 3 x(n(x)-Yl t*rNr) Here pj (j — 1,2,3) are the chemical potentials associated with the lepton, baryon 

and electric charge conservation, so that p n = fj,2, A*p = Pi + (^3,p e = A*i + M3 an d = Pi- The thermodynamical 
potential can be written as f2 = —f3~ 1 lnZ . Usually the eigenvalues of TL contain the contribution from neutrons, 
protons, electrons and some meson species, and the densities are N r = —d£l/dp r , where r — n,p,e... We name 
n = J2 r tt r , M. — ^2 ir M. r the total thermodynamical potential and magnetization, respectively. 

A standard procedure is to work in the mean field approximation in which the meson fields a, p, lu are taken as 
constant, as done in Refs.|23|. [24|. through which the mass spectrum of baryons is corrected and strong repulsive 
interactions between them is found. However, for simplicity we will keep the spectra in the tree approximation to 
obtain the one-loop approximation for f2, and neglect the statistical contribution from meson terms in Q as compared 
with those of fermions (for them niif3 ~ 10 3 — 10 4 ) except for fields B < B cp = m? p /e ss 10 20 G, since the contribution 
of the p vector meson condensate to M. becomes relevant, and in analogy with H /± -sQ, leads to a self-consistent 

spontaneous magnetization B = inAi = 2ireN p ^Jm 2 — eB, where N p is the condensate density. However, for such 

fields the magnetic pressure: J AidB — M.B = —B 2 /8ir overwhelms the kinetic pressure term P (of order 10 36 
dynes/cm 2 ) leading to P± < 0, and the star is definitely unstable: it collapses. This mechanism is valid for other 
quasi-particle vector boson condensates, as di-quarks, which may be formed in the medium even for smaller values of 
the field B 



THERMODYNAMICAL POTENTIAL OF A NEUTRON GAS 

The Green functions method [2^| leads to a general expression for the relativistic thermodynamical potential. At the 
one-loop level, where no radiative corrections are considered, it is a generalization of the usual nonrelativistic formula 
because of the fact that antiparticles must also be included. Particles and antiparticles contribute with chemical 
potentials of opposite sign, leading to sums or integrals over the quantum numbers involved, of terms containing the 
product of the logarithms of (e _ ( M " ±B ™) + 1). In this case, when no external fields are present, a divergent term 
accounting for the vacuum energy appears which is usually subtracted j^. In presence of an external field, however, 
a term accoun ting for the vacuum contribution, must also be included: The Euler-Heisenberg energy of vacuum in an 
external field [ill Hil l30t l3l| . One can obtain an expression for the thermodynamic potential of the neutron gas in 
the one-loop approximation as Q n = Q sn + fiy n , with 

1 f°° 

^ sn = ~T~27i z2 / P±dp±dp 3 ln [f + (p n , (3)f~ (p„, P)] , (13) 



where ^(p^P) - (1 + e - ^"^")^) accounts, respectively, for the contribution of particles and antiparticles. The 
expression for the vacuum term reads thus 

1 f°° 

ttvn = ^2 I P±dp ± dp 3 E n (14) 



'/= 



which is divergent. In the Appendix we will show how to regularize this expression, and how to obtain the analog to 
the Euler-Heisenberg energy of vacuum due to the neutrons contribution. 

After integrating by parts in (|13fl its evaluation becomes easier. The Fermi distributions, which arise by differ- 
entiating / ± with respect to pz, are = 1/(1 + e^ En ^ fJ ' n '^). In the degenerate case this expression reduces in 
n~ = 8(p — E n ) and n + = 0, since in that case only particles contribute to f2. The resulting expression splits itself 
in two terms where the integrals are bounded by the Fermi surfaces p — E n (jj = ±1) = 0. These surfaces have axial 
symmetry, and thus, the Fermi momentum is not a definite number, given only in terms of p n and m n , but on the 
opposite, it has infinite values. Thence, we have 

—A E r 'p±dp ± ^e(p-E n (r 1 )) . (15) 

77=1,-1 *' Q n 



The 9 function bound these integrals in the intervals —p^F P3 P3F, where p 3 p = y p 2 — {\J p\ + rrv^ + r/y) and 
< p± < y/(jt — r]y) 2 — m 2 . After some transformations, it yields 
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-no E 

77=1,-1 

(1 + W) 3 (3 



12 

- m) 



(1 + r}y)(5r)y - 3)xf 



24 



24 



777/2; 

~6~ 



(16) 



where a; = fx n /m n , (m n (x — 1) is the usual Fermi energy), and y = qB/m. We define the functions / = f(x,r/y) = 
\J x 2 — (1 + ?7j/) 2 , s = s(x,r/y) = (tt/2 — sin -1 (l + rjy)/x), L = L(x,r/y) = In (a; + f(x,rjy))/(l + r/y). The functions 
f =< Pf > /m n , where < >= Wfj^— (m n + rjqB) 2 and we name them the average Fermi momenta for rj = ±1. 
We see that m n — ?/<7-B behaves formally as a two-valued magnetic mass. 
In the zero field limit one gets Pq = —Vl n (y = 0), where 



tt sn (y = 0) 



-12~ — + 8 Lo 



where /o = \/ x 2 — 1 is the relative Fermi momentum pp/m, and Lo = ln(a; 
Appendix) has an Euler-Heisenbcrg-likc form as 



(17) 



/o). The neutron vacuum term (see 



r2y n = 



1 

4^2 



dyy e 



-3„-(m„+?"-B")a 



cosh(gi?r7iy) - 1 - (qBmy) 2 /2\] 



(18) 



+ ^-5- / dyy" 2 / du;e- [( "' +m " )2+ « 2 - B2ly [sinh(2 (7 B(ii; + m„)y) - (2gB(u; + m„)y) + (2gB(w + m„))/) 3 /3!] 
2?r Jo Jo 



It can be shown (see Appendix) that the more significant term in i|18|) is the first one, which for fields of order 10 17 
G leads to ttvn ~ 10 30 erg-cm -3 and is negligible small as compared with fl sn up to B ~ 10 18 G. Thus, we neglect 
it in a first approximation in what follows. We must point out, however, that since neutrons have a quark structure, 
a more fundamental quantity would be the vacuum quark contribution, whose order of magnitude is expected to be 
near f2y„. Apart from this note, it should be emphasized that the role of vacuum cannot be ignored for fields greater 
than 10 18 G. 

From N n = dfl/d^ n one gets 



N n = N 

77=1,-1 

In the limit B — Eq. I|19H reproduces the usual density of a relativistic Fermi gas at zero temperature, N n — Nof§ /3 
Having an equation relating the chemical potentials, and demanding conservation of both baryonic number N n + 
N p = Nb and electric charge N p + N e — 0, in principle one may solve exactly the problem in terms of the external 
field as a parameter. Nonetheless, we shall focus our discussion on the properties of the equation of state. Note in 
passing that our expressions for the spectra and densities of neutrons and protons are similar to those of Ref. of 
a neutron gas in a magnetic field, but we get different equations of state. 
Finally, for the magnetization, given as M. n = —dfl n /dB, we have 



t 
3 



vvO- + vy)f) vyx 



(19) 



-Mo V 

77=1,-1 

(l+ W ) 2 (l-?7t//2) r , x 



(1 - 2vy)xf 
6 



(20) 



L 



6 



where A^ = m 3 /47r 2 - 2.04 x 10 39 , Q a = N m n - 3.0 x 10 36 , and M = N a q ~ 2.92 x 10 16 and one can write 
M n = M + (rj = -1) - M+(rj = +1), and obviously, M n > 0. 

We confirmed by explicit calculation that M is a nonlinear function of B and, in this sense, the magnetic response 
is ferromagnetic. A fully ferromagnetic response demands to include also the spin-spin coupling contribution. We 
discuss briefly this point below. 
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To see why the magnetization is always positive for the neutron gas note that the magnetic susceptibility x = 
dA4 n /dB can be easily obtained as 



which for x > 1, y < 1, and / real and positive, it is x > 0. This means that M n is an increasing function of B (or y) 
under such conditions. As M n (y = 0) = and M n (y = 1) = 2Mv{l + vr) > 0, this means that M n > in the region 
between these two points, which is the one of interest for us (the region I discussed below and showed in Fig. 2). 

The fact that we are summing over the magnetic moments oriented parallel (?y = —1) and antiparallcl (77 = +1) to 
B is similar to the well known Pauli paramagnetism in nonrelativistic quantum statistics. We may consider each term 
i] = ±1 as representing a phase of the system. One can draw in the x,y plane, for both x > 0, y > 0, two regions 
limited by the lines x — y = ±1 and x + y = 1 (see Figure 1). The region I is bounded above by the line x = y + 1, 
on which the contribution from f(rj = +1), and all other terms containing r\ = +1 vanish: all magnetic moments are 
aligned parallel to B. Below, such region is limited by y = 0. For points in the region x > y + 1 and y > both terms 
containing r\ = — 1, +1 are nonzero and real. Thus, the quantities tt, N n , and M. are real and although most neutrons 
have their magnetic moments along B, there are some amount of them having their magnetic moments antiparallcl 
to B. In this region, near the line x = y + 1, is located the curve P±_ = 0, which we will discuss below. 

In the region II limited by x > 1 — y, y + 1 > x > y — 1 only the term 77 = — 1 is real and the term 77 = +1 becomes 
purely imaginary, and all neutrons would have their magnetic moments parallel to B. The term 77 = —1 contribute 
to Q, N n , and M. For x < y — 1, there are no physical solutions. For y = 0, the magnetic field, and in consequence, 
the magnetization M vanish. For y < (magnetic field along the negative y axis) we have the reverse situation to 
the one described above (See Figure 1). 

With regards to the background of electrons and protons, we remind that the case of the electron gas was discussed 
m (0D- We expect that for magnetic fields of order B ce — m 2 /e ~ 10 13 G and densities around 10 30 cm -3 all 
electrons are in the Landau ground state, and the system show the instability which arises from the vanishing of the 
transverse pressure. For densities and magnetic fields above these critical values, the stability of the electron gas is 
doubtful. Recently Ghosh, Mandal and Chakrabarty |j| has pointed out strong arguments against the equilibrium 
of an electron gas under f3 decay in strongly magnetized neutron stars. The situation may be modified if radiative 
corrections are taken into account: if B ce = m 2 e /e ~ 10 13 G is the usual QED critical field, for larger fields of 
order B' ce ~ 4ir /aB ce ~ 10 16 G the contribution of the electron anomalous magnetic moment becomes significant 
and the problem cannot be satisfactory treated at the tree level. However, starting from the results of uj,|£j|, one 
concludes that the electron gas is hardly in equilibrium for fields beyond B ce . One possibility is the bosonization 
of the electron system, as has been recently suggested j3j| • This may be accomplished through the increase of the 
spin-spin interaction, leading to the formation of parallel spin electron pairs, equivalent to charged vector bosons. 

If we include both the normal and the anomalous magnetic moments for protons, one can give a formula for their 
spectrum in the external field B as|l4j: 



q p = 2.79M n . For neutrons, the critical field at which the coupling energy of its magnetic moment equals the rest 
energy is B cn = 1.57 x 10 20 G. For protons it is B cp = 2.29 x 10 20 G. By defining x p — jJL p /m p , y p = q p /m p , b — 2e/m 2 , 



then y p = 2.79e/2m 2 . We name also g = g(x p , B, n) = yj x p — h(B, n) 2 and h = h(B, n) — (\JbBn + 1 + -qy p B). Thus 

for the proton thermodynamical potential we get 




(21) 




(22) 





n ±77 



(23) 



and for its density 




n ±7] 



(24) 



10 



FIG. 1: This plot shows the regions in the x,y plane where the neutron magnetic moments are oriented parallel or antiparallel 
to the magnetic field B. Special attention should be given to regions I and II, where the solutions discussed in the text are 
valid. 



while the magnetization is given by 



emz 



n ±7/ 

+ (bn/2y/bBn + 1)) x In (x p + g) /h\ , 

where the coefficients of these formulae are N Qp = em p B/2ir 2 ~ 4.06 x 10 19 5, n 0p = N Qp m p B ~ 6.1 x 10 16 5, and 
■Mop = No p m p = Hq/B. The maximum occupied Landau quantum number n may be given as n max = (x p — r]y p B) 2 — 
1/bB. For B <C B cp , so that y p B <C 1, and x p > 1, one can take approximately n max ~ (x 2 — l)/bB, and for fields 
large enough n max = 0. We expect that from the equation fi n = ji p + fj, e , then x p ~ x n , the previous expression for 
the proton density N p decreases with increasing B, favoring the inverse beta decay. For fields B ~ m p /q p and x p 3> 1, 
n max > 1; and thus large Landau numbers are again occupied. However, for x n ,x p > 1, being both quantities of the 
same order of magnitude, from the comparison of Qq p , No p ,M.Qp with Qq, No, Mo we conclude that for fields below 
10 19 G, the dominant longitudinal pressure, density and magnetization comes from the neutron gas. 



CONDITION FOR ZERO TRANSVERSE PRESSURE AND COLLAPSE 



In the electron gas case |l| , the vanishing of the transverse momentum can be guessed from the spectrum when all 
the system is confined to the Landau ground state. The spectrum corresponds to a purely one-dimensional system 
moving along the external field, and the transverse Fermi momentum is zero. In a similar wayin the neutron gas case 
we observe that the threshold of zero transverse pressure P± = can be figured out from the spectrum described by 
Eq. (|12fl . since the contribution from the T] = — 1 term is dominant (observe that the term with rj = +1 contributes 
with a negative term to A4). We shall consider on the Fermi surface for rj = — 1 the quantity 



PF±ef f = Mn - PF3 ~ m n = [V Pf± + m l ~ V B ) ~ m n (26) 

which we name the effective squared Fermi transverse momentum. If B <C 10 20 G, then q 2 B 2 <C 2eBm n , the vanishing 
of PF±eff is guaranteed if 



p F ± - ^2qBm n . (27) 

The resulting Fermi surface would be equivalent to that for one-dimensional motion, parallel to B, i. e., for particles 
having energy E n ~ \/v\ + TO n: an d in consequence the transverse momentum (and pressure) vanishes. Notice from 
(12711 that for pF±/rrin ~ 10~ 1,5 one has qB/m n = y ~ 10 -3 , which means fields of order 10 17 G. A more accurate 
result is obtained, however, from the condition: 7j_ = 0. 

In Figure 2 we have drawn the equation P± n = — ft — BAi = in terms of the variables x, y. We observe that there 
is a continuous range of values of x, starting from x = 1.005 to x = 1.125 and from y = 0.001 to y = 0.1, for which 
the collapse takes place. The latter values of y means fields in the interval B ~ 10 17 to 10 19 G. To these ranges of x, y 
corresponds a continuum range of densities, from 10 _2 A^o (10 12 g/cm -3 ) onwards. The transverse compression of the 
whole mass of the star due to flux conservation leads to an increase of B and the mechanism of collapse is enhanced. 



The spin-spin coupling contribution 



Although we used previously the term ferromagnetic to denote the non-linear response of the medium to the 
external field, what we have considered actually in our previous calculations is the occurrence of relativistic Pauli's 
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FIG. 2: The instability condition: P± = in terms of the variables x, y. A nascent strongly magnetized neutron star having 
a configuration such that its dynamical stage would be represented by a point above the central curve in this plot would be 
unstable to transversal collapse, since P± < there. 



FIG. 3: The curve P± = in terms of the neutron star mean density N and the surface magnetic field B. 



paramagnetism. We have not fully considered the spin-spin coupling, which would lead more definitely to Heisenberg 
ferromagnetism, and is physically reasonable to expect its appearance in nuclear matter for densities high enough. 
This problem has been studied in [34| via the interaction through axial vector and tensor exchange channels. These 
authors show that if the interactions are strong enough and differ in sign, the system loses the spherical symmetry 
due to a mechanism independent of the one discussed in the present paper. In the case of a quark liquid |35l | the 
problem has been studied under one gluon exchange interaction, and conditions for ferromagnetism may arise. 

Although the problem requires further research, one expects that if Heisenberg ferromagnetism is to appear, it 
would increase the magnetization to M s — kM., where k is the internal field parameter. If k » 1 our previous 
estimates for the vanishing of the transverse pressure might be largely exceeded, with the arising of a new spontaneous 
magnetization M s ~ kM. = B/4tt. This would mean that the magnetic field B could be kept self-consistently and 
our previous calculations would be a lower bound of the fully ferromagnetic case. The vanishing of P± is expected 
to occur surely at values of B smaller than those depicted in Figures 2 and 3. As a rough estimation, if we assume 
the exchange interaction J among neutrons of order of few hundreds of MeV, and the number of nearest neighbors as 
z ~ 10, by dividing their product by the dipole interaction energy, say, for the core of the star where N = N , we get 
k ~ 10 4 . This means B ~ 4irA4 s ~ 10 20 G. For such extremely large fields the magnetic coupling of quarks with B 
would become of the order of their binding energy through the color field producing a deconfinment phase transition 
leading to a quark (q)-star, a pressure-induced transition to uds-quark matter via ud-quark condensates, as discussed 
in Refs. |43l l44j . But fields of that order lead surely to the collapse of the star and even to the instability of vacuum. 



MAGNETAR FORMATION AND STABILITY 



Next we briefly review the basic ideas supporting the theory of magnetars and then show why these hypothetical 
objects cannot survive after reaching the claimed super strong magnetic fields. We then present prospectives for 
a hybrid or strange star to appear as a remnant of the quantum magnetic collapse of a NS. According to Duncan 
and ThompsonQ, NSs with very high dipole surface magnetic field strength, B ~ [10 14 — 10 15 ] G, may form when 
(classical) conditions for a helical dynamo action are efficiently met during the seconds following the core-collapse and 
bounce in a supernova (SN) explosion 3] . A newly-born NS may undergo vigorous convection during the first 30 s 
following its formation 39] . If the NS spins (differentially) sufficiently fast (P ~ 1 ms) the conditions are created for 
the a — Q dynamo action to be built, which may survive depletion due to turbulent diffusion. Collapse theory, on the 
other hand, shows that some pre-supernova stellar cores could adquire enough spin so as to rotate near their Keplerian 
equatorial velocity, the break-up spin frequency: Qk ^Illlf GnM/R 3 ) 1 / 2 , which implies a period Pk ~ 0.6 s, after 
core bounce. Under these conditions, fields as large as |3j,|42j 



B~\tf?(JL\ G, (28) 
\ 1ms J 

may be generated as long as the differential rotation is dragged out by the growing magnetic stresses. For this 
process to efficiently operate the ratio between the spin rate (P) and the convection overturn time scale (r con ), the 
Rossby number (Ro), should be Rq < 1. Duncan and Thompson warned that JJo » 1 should induce less effective 
mean-dynamos 3J). In this case, an ordinary dipole B D ~ [10 12 - 10 13 ]G may be built by incoherent superposition 
of small dipoles of characteristic size A ~ [| — l]km, and a saturation strength B sat — (inpY^X/Tcon ~ 10 16 G may 
be reached at the surface during this early evolution of the nascent neutron star. At such fields, the huge rotational 
energy of a NS spinning at: cjns > 1kHz, is leaked out via magnetic braking and an enormous energy is injected into 
the SN remnant. This process may explain the power emitted by a plerion. 
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As shown above, at the end of the SN core collapse we are left with a rapidly rotating NS endowed with with an 
extremely strong magnetic field (ESMF) strength and a large matter density p ~ [10 14 — 10 12 ] g cm~ 3 . As illustrated 
in Figure 1, those are the conditions for the quantum instability to start to dominate the dynamics of the young pulsar. 
At this stage, the magnetic pressure inwards may overpass the star energy density at its equator and the collapse 
becomes unavoidable. As the collapse proceeds, higher and higher densities are reached till the point the supranuclear 
density may reverse the direction of implosion. A hybrid or strange star (SS) may have formed. We explore next 
this plausible outcome, among the other possibilities quoted above. From that moment, the sound wave generated 
at the core bounce builds itself into a shock wave traveling through the star at Vsw ~ c/V3 km s" 1 . Although 
the ESMF strength could be quite large as long as the collapse advances, the huge kinetic energy, E ~ io 51 ~ 52 erg, 
the mean energy obtained in calculations of energy release in neutron star phase transitions to strange {twin) stars 
pjil I37I l3S| and some prompt shock supernovae {E p _ s ~ 10 51 erg) |3^|, carried away by the shock wave drives a kind 
of supernova explosion inasmuch as in the quark nova model and similar scenarios 38]. Such a huge ram pressure 
may counterbalance the magnetic pressure, and even surpass it, i. e., 

B 2 { R \ 6 

PejectVl w >- — , (29) 

8717x0 \rAJ 

at a location from the star center equivalent to the Alfven radius of the magnetar 



Ta 



/2tt 2n 


f 


"B 4 i? 12 " 






MM 2 . 



1/7 

- 80 km. (30) 



This radius is quite large, about 7 times the NS radius (see Table 1). Here M defines the accretion rate of the 
free-falling overlaying material making up the NS crust, which is left out when the transition occurs (see definitions 
and further details in Ref.[3^], and references therein). Therefore, it is quite legitimated to expect that most of the 
magnetic energy stored inside the magnetosphere to be drained out of the Alfven radius. Notice in addition that the 
strange star radius scales as: i?gs = -Rns(*~) ; where the densities ratio reads: ~ 0.1 — 0.2. Other relations 
between both the stars can be obtained by using conservation laws or appropriate scalings. 

Then the ESMF lines are pushed out and finally broken, in a process inverse to the standard accretion one, from 
onwards into the SN remnant surroundings, as a violent explosion that dissipates a large part of the magnetic flux (<fr ~ 
B 2 r 2 A ) and energy trapped inside the magnetar magnetosphere |2g. Energy from the magnetic field can be dissipated 
via vacuum polarization and electron- pos itron pairs creation, as well as acceleration of charged matter flowing away 
(synchroton and curvature losses |4fl l4lj with the explosion and material trapped in the star magnetosphere and 
nearby the Alfven radius. This is analogous to the mechanism operating during a solar flare or a coronal mass- 
ejection, where the very high B in a given Sun-spot is drastically diminished after flaring for a short period of time 
(see also Ref. ^3)- ^ n * ne Sun s P°ts outbursts and coronal mass ejections launch into space part of the Solar wind of 
charged and neutral particles passing by the Earth^^,^3- I n t ne case °f an imploding NS, the phase of open magnetic 
field lines over which the strange star is acting as a propeller lasts for about AT prop ~ E spin / L prop ~ [10 2 — 10 3 ] s, 
with L prop r~j 2Mc 2 the propeller luminosity, and E sp i n ~ IcJ$$ the star rotational energy. Thence, the large amount 
of matter ejected from the strange star at such large velocities and the pairs created, in the vacuum breakdown, drains 
out the dipole field of the remnant below the quantum electrodynamic limit of B ce ~ 4.4 x 10 13 G 28]. 

To give an insight into this piece of the physics of the problem, notice that once the propeller phase is over and no 
more luminosity is coming from that mechanism, the magnetic field lines can recombine again if the energy released 
in this new stable phase is essentially the strange star rotational dipole luminosity (as the one from a millisecond 
pulsar), which then becomes the star dominant mechanism of energy emission, that is 



4 g ssgs^s =lLJss cJ ss . (31) 

For the parameters shown in Table 1 and the observed luminosity from fast rotating pulsars, this relation implies 
a new equilibrium magnetic field B ~ 10 12 ~ 13 G, which is below the quantum electrodynamics threshold. 

Since all the differential rotation has been dragged up to build up the former ESMF, then nothing else remains to 
make the magnetic field to grow to its pre-collapse value. Thence no such ultra high B should reappear. We may be 
left with a sub-millisecond strange star^^] or a hybrid star0] with "canonical" field strength, but no any magnetar. 
We note in passing that the above theoretical result is attained on the standardized assumption that the structure of 
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the magnetic field of the pulsar is dipolar. This premise is underlying to the claim by Kouveliotou et al.[4{j that a 
magnetar had been discovered in the soft gamma-ray repeater source SGR 1806-20. Notwithstanding, for other NS 
(multipolar or uniform) field configurations we do expect the overall behavior here discussed to persist, since once 
the spin-spin coupling is taken into account the unavoidable consequence is the appearance of a ferromagnetic (axial) 
configuration or structure which is dipolar in nature, and therefore the theory propound here still holds, because the 
magnetic field could then be amplified by a factor 10 4 ), putting the nascent pulsar above the threshold for stability, 
and the collapse ensues. Rephrasing this, one can think of this theory as a field-configuration independent constraint 
on initial NS magnetic field strengths. 

In looking for other contexts far from those involving compact remnant stars, we noticed that recent "Tabletop 
Astrophysics" experiments performed by C. Wieman et al.[42j have succeeded in refining sophisticated techniques 
to switch atoms in a Bose-Einstein condensate (BEC) |48| from states of implosion to states of re-expansion (or 
explosion). In certain cases they observed that some collapses appear rather similar to microscopic supernovae. The 
initial implosion is followed by an explosion in which atoms are ejected as a hollow ball or in narrow jets, like in 
the collapse and rebound of a exploding star that forms characteristic expanding balls or streams of outflowing gas. 
The Wieman team named the phenomenon "Bosenovae" because of the similarities with typical supernovae. In fact, 
in both phenomena some material is left after the implosion as a compact object. Because of the similarity of the 
physics of Bose-Einstein condensates with the one we presented above, we are confident that this our theory can also be 
adapted to explain these impressive results by Wieman et al.|4^; Ketterle and Anglin 26] . The main feature of these 
oscillating Bose-Einstein condensates, said to be a scale-down version of either a neutron star or a white dwarf [4^, is 
that for some critical fields what appears is an attractive force between the atoms and the condensate implodes and 
rebounds driven by some sort of internal negative pressure. We claim in this paper (a detailed description of BECs 
phenomenology is to be given in a work in preparation^^) that such a negative pressure could be explained in the 
context of the theory introduced in this paper, since a precise relationship between gas density and magnetic field 
strength in the BEC is settled out by switching the atoms between attractive and non-attractive states. Thence the 
negative pressure acts as the equivalent of an attractive force among the atoms directed towards the magnetic field 
axis, leading to the BEC implosion. 



CONCLUSIONS 



We conclude by claiming that if a degenerate neutron gas is under the action of a super strong magnetic field 
B C e P$ B < -B C n, for values of the density typical of NS matter its transverse pressure vanishes, the outcome being 
a transverse collapse. This phenomenon establishes a strong bound on the magnetic field strength expected to be 
found in any stable neutron star pulsar, regardless of its initial field configuration, and suggests a possible endpoint 
in the early evolution of highly magnetized neutron stars. They could likely be a mixed phase of nucleons and 
a 7r^~ , 7T°, K ± ,K°,K° ,a,pr^,u} meson condensate, a hybrid or strange star, or a distorted black hole but no any 
magnetar at all. We point out, nonetheless, that if by any mechanism a strange star could be formed directly in a 
supernova cxlposion (which is uncertain, but not ruled out) and if vigorous dynamo action operates in the strange 
quark matter bulk, "pulsars" with fields higher than B ce could still be formed. In other words, if any magnetars exist, 
they cannot be neutron stars. 



APPENDIX 



We use the integral representation 



a 1 ' 2 = 7T" 1 / 2 r %{e~ ax * - 1) = 7T- 1 / 2 / dyy^' 2 {e- a y - 1). (32) 



x 



We regularize the divergent term dependent on a in l|32(l . by introducing a small quantity e as the lower limit in 
the integral, and neglecting the term independent of a 



,1/2 



r-V2 



dyy 



-3/2 e -ay_ 



(33) 
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By taking a(e, 77) = p 2 + ( \J p\ + m 2 + rjqB) 2 and substituting in (|15(l and performing the Gaussian integral on p%, 
one obtains 



n Vn (e) = V / p±dp± / dyy-5/2 e -fV^i+^+wB)y (34) 



By substituting 2 = -^/p j_ + ?77^ + rjqB one is left with the expression 



1 r 

n vn(e) = ^ / d 2/2/ - 3 e-( m "+' ?2B2 ^cosh( (Z Bm2;) 

?7=±1 47r Jm n +7jqB 



By introducing the new variable w = z — m„ — ?7(?i?, the second integral in l|35[) becomes 

USE [°° AMI -a /■°° p -[( 10 + m „) 2 + 9 2 B 2 ] a 



4tt 2 



dyiT 2 / e -l^+ m -J +9 " I* sinh(2 9 B(u; + m„)j/). (36) 



After subtracting to cosh(qBmy) and sinh[2q_B(w + m n )y] the first two terms in their series expansion, one can 
take e — > and obtain the finite expression (|18fl . This process is equivalent to the subtraction of divergent terms, one 
of which is proportional to B 2 , and absorbs the classical field energy term B 2 /8tt. 

It is not difficult to check that for fields B <C 10 20 G, the first term in 118(1 is the dominant one. Its first contribution 
after the series expansion of cosh(qBmy) is q 3 B 3 m^/2TT 2 (m 2 l + q 2 B 2 ). For fields of order 10 17 G such a term is of 
order 10 30 ergs/cm 3 , much smaller than Sl sn . But for fields near 10 20 G its contribution is comparable to that of VL sn . 
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TABLE I: Parameters used in modeling the phase transition of a neutron to a strange star, as discussed in the text. 
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